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Abstract. - We investigate how the topology of small-world networks affects the dynamics 
of the voter model for opinion formation. We show that, contrary to what occurs on regular 
topologies with local interactions, the voter model on small- world networks does not display the 
emergence of complete order in the thermodynamic limit. The system settles in a stationary 
state with coexisting opinions whose lifetime diverges with the system size. Hence the nontrivial 
connectivity pattern leads to the counterintuitive conclusion that long-range connections inhibit 
the ordering process. However, for networks of finite size, for which full uniformity is reached, 
the ordering process takes a time shorter than on a regular lattice of the same size. 



In the last decade, social sciences have started to deal with large scale modeling of a variety 
of spreading and ordering phenomena that involve cooperative behavior 1 . In this context, 
classical models developed in statistical physics to study the onset of order in matter j2| have 
turned out to be useful for the investigation of the principles at the basis of social ordering. 
For instance, the Ising and voter models and their variations are prototypical models for a 
wide class of social interaction phenomena |2101|51|Sl|71|Sl|Uj . The voter model is possibly the 
minimal model for opinion spreading and the study of the onset of consensus. It is usually 
defined on a regular lattice of dimension d. Each site is characterized by a discrete variable 
s that may assume two values (s = ±1) representing two opposite opinions; for instance 
the electoral choice in favor of two different candidates. Starting from a disordered initial 
configuration, the model follows a simple dynamical evolution in which at each time step one 
site is selected at random and made equal to one of its nearest neighbors (chosen at random 
on its turn). This dynamics mimics the homogcnization of opinions through the confrontation 
of peers and leads to the formation and coarsening of ordered regions where individuals share 
the same opinion. In d = 1 and 2 the model eventually converges to an ordered state with all 
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variables having the same value |1(J| . This state is absorbing since the system cannot escape 
from it once it is reached 

While the use of regular lattices to model the interaction between elementary objects is 
well justified for most physical situations, such an assumption is not obvious in the context 
of social sciences. Many social systems indeed show interaction patterns that find a better 
characterization as complex networks with distinctive connectivity properties |12lll!'{| . In a 
graph representation, where nodes identify the individuals and links their direct interactions, 
many social and natural networks exhibit peculiar topological properties related to the pres- 
ence of highly connected individuals and long range connections. Among these features, the 
most well documented is the small diameter of social networks, i. e. each individual can reach 
any other one passing through a very small number of intermediate nodes. In addition, social 
interactions favor the connection between common acquaintances leading to the presence of 
high clustering among nodes. This is quantitatively expressed by a high probability that if 
two nodes share a neighbor they are directly connected on their turn. This last property, 
along with the small diameter, define the so-called small- world behavior 14[ll5j . The use of 
small-world like topologies in models of opinion spreading is then a logical step in the direc- 
tion of a more realistic approach to the phenomenon. The interaction patterns involved in 
the process of opinion formation are very likely similar to those of social networks such as the 
web of sexual contacts JJ or scientific collaborations |17| . 

In order to investigate how complex connectivity patterns might influence opinion forma- 
tion, we study the effect of the small-world topology on the evolution of the voter model. The 
prototypical network possessing the small-world character is the Watts-Strogatz (WS) model 
that has been extensively studied in several contexts |14[ll5l[TS| . In particular, we have con- 
sidered the WS network as defined in Ref. Starting from a one-dimensional lattice of N 
sites with periodic boundary conditions and each node connected with 2k nearest neighbors, a 
stochastic rewiring is introduced. Nodes are visited one by one sequentially and each of the k 
links connecting the node to its nearest neighbors in the clockwise sense is rewired with prob- 
ability p to a randomly chosen node. As p is increased, the WS network interpolates between 
a one-dimensional lattice {p — 0), with only geographical neighbors in contact, and a random 
graph {p = 1), where short and long range connections are equally likely. The small- world 
behavior (small diameter, high clustering) is exhibited for values of the rewiring probability 
p such that \/{kN) <C p ^ 1. The transition between the one-dimensional topology and the 
small- world one occurring for p « l/{kN) is governed by the value of ^, the average distance 
between nodes connected with shortcuts. ^ is the only nontrivial length in the network and 
can be simply shown to scale as l/{kp) [19II20| . If the network size N is much smaller than ^, 
that is p <C l/{kN) the system does not have long-range connections and is a one-dimensional 
lattice. For N ^ ^ many shortcuts are present and originate the small-world behavior. 

For the study of the ordering of the voter model, the natural quantity of interest is the 
fraction ua of active bonds, i. e. the density of links connecting sites with opposite values 
of s. These are the links where the dynamics takes place. Such a quantity vanishes if the 
system orders entirely, while it remains finite if domains coexist, and I/ha is a measure of 
the average size of domains. In the case p = (one-dimensional system), due to the diffusive 
motion of interfaces between domains, the fraction of active bonds decays with respect to time 
t as UA 123 up to a crossover time ^ N"^, after which ua exhibits an exponential 

relaxation to the absorbing state = (see Fig. Such a fast decay is the effect of the 
finite size of the system. In the thermodynamic limit one recovers a pure power-law relaxation 
to the absorbing state. 

We consider now a Watts-Strogatz network with k = 2, p = 0.05 and values of N such 
that the system is safely inside the small- world regime {N > 200). The behavior changes 
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Fig. 1 - Log-log plot of the fraction ua of active bonds between nodes with different opinion. Values 
are averaged over 1000 runs. Time is measured in Monte Carlo steps per site. Empty symbols are for 
the one-dimensional case (p = 0). Filled symbols are for rewiring probability p = 0.05. Data are for 
iV = 200 (circles), iV = 400 (squares), iV = 800 (diamonds), iV = 1600 (triangles up) and iV = 3200 
(triangles left). 



dramatically (Fig. 0: After a transient, the plot of tia exhibits a plateau, indicating that 
domains remain on average of constant size. This regime is ended by an exponential approach 
to the absorbing state ua = 0. The duration of the plateau grows with N , but (as can 
be seen from Fig. ^) the time to reach complete ordering is smaller for the WS network 
than for a regular lattice with the same number of nodes. This is in agreement with the 
naive expectation that long-range connections should speed up the homogenization process. 
However, this occurs in a highly nontrivial way: During most of the evolution, ua is higher on 
the small-world network than on a regular lattice, i. e. the small-world network is, for a long 
time interval, more disordered, and orders rapidly only at the very end. We have checked that 
this phenomenology is not an artifact of the rewiring procedure to build the Watts-Strogatz 
network: analogous results are obtained when the small-world topology is produced by adding 
random links to a one-dimensional lattice. 

The nature of the exponential approach to the absorbing state is related to a standard 
finite-size effect in the presence of absorbing states. Any finite system settles in a stationary 
state with constant activity until it hits the absorbing state because of a large spontaneous 
fluctuation [TTl- The survival probability Ps{t) that the system is still in an active state after 
time t decays exponentially, Ps{t) ~ exp(— t/r). Here t is the average lifetime in the active 
state and is found to increase with the system size as t ^ (Fig. |21 inset). This implies that 
in the thermodynamic limit (A'^ cxd), the system remains indefinitely in the stationary state, 
with everlasting activity, i. e. incomplete ordering. Hence, the voter model in the small-world 
regime behaves as its mean field version (euclidean lattice with d = oo) that does not reach an 
ordered state |2H . This finding is quite interesting: in the thermodynamic limit, the presence 
of long range connections does not make the ordering process easier, rather it inhibits it. The 
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Fig. 2 - Main: The survival probability Ps{t) as a function of time for a network with p = 0.05. Data 
are for A'' = 200 (circles), A'' — 400 (squares), A'' = 800 (diamonds) and A'' — 1600 (triangles) and are 
computed averaging over 1000 runs. Inset: The average lifetime in the active state r as a function of 
the system size A'. Circles are numerical values obtained by measuring r as the inverse slope of the 
curves shown in the main part of the figure. The solid line is the best fit to the expression r cx A'^'' 
giving 7 = 1.02 ± 0.02. 



small-viforld topology of the network represents a barrier against convergence to order. 

In order to understand better the origin of this incomplete ordering, we study networks 
with N = 10^ nodes, for which the lifetime r is much larger than the time scales of interest 
and may therefore be considered in practice as infinite. In this case, as long as p > lO^'"^ the 
network is in the small world-regime and ua tends to a finite stationary value (see Fig. |3J)- 
The network settles in a dynamically active regime in which complete consensus does not 
emerge. The value of ua in the stationary state depends on p. In the one-dimensional lattice, 
the ordering process takes place via free diffusion of domain boundaries (active bonds) and 
their annihilation upon encounter. In the small-world networks, shortcuts are an obstacle 
for free diffusion of active bonds and tend to pin domain boundaries. More in detail, the 
role played by sites connected with shortcuts underlies the following simple scaling analysis, 
valid for t ^ 1 and small p. For short times domains form and start to coarsen: I/ua, their 
average size, is much smaller than the average shortcut distance ^ = 2/{kp). In this regime 
the evolution is practically equal to the one-dimensional voter model, with ua decreasing as 
When the length of domains reaches the size of ^, the only other length in the problem, 
the behavior changes. The cross-over to the stationary state occurs therefore when I/ua ^ ^, 
defining a diverging crossover time t* p^^- Since the crossover takes place at ua ^ p we 
obtain that 

nAit,p)=pgit/p-^), (1) 

where the scaling function g{x) approaches a constant value for x ^ 1. Equation is well 
obeyed by numerical results, that show a good data collapse on the predicted behavior (Fig. 13 
inset). 
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Fig. 3 - Main: The fraction ua of active links as a function of time for a large network with A*" — 10^ 
nodes. Data are for p = 0.0002 (circles), p — 0.0005 (squares), p — 0.001 (diamonds), p = 0.002 
(triangles up) and p — 0.005 (triangles left) Inset: the same data, divided by p and plotted as a 
function of t/p"^ in order to show the validity of Eq. iQ. 

In summary we have shown that complex topological properties of small-world networks 
strongly affect the behavior of the voter model, leading to incomplete ordering in the thermo- 
dynamical limit and to counterintuitive phenomena for small systems. We believe that the 
behavior of more general models for social influence [22] is similarly modified by the topol- 
ogy of the interaction network. It would be interesting to test on these models the effect of 
more heterogeneous topologies, such as scale-free networks, which are known to alter several 
dynamical processes occurring on them PSIEHES] . 

During the completion of this work, we have become aware of some recent work |2fi| on 
the ordering process of the Ising model on the WS network, presenting conclusions somewhat 
similar to ours. 

* * * 
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REFERENCES 

[1] AxELROD R., The complexity of cooperation (Princeton University Press, Princeton) 1997 
[2] Bray A. J., Adv. Phys., 43 (1994) 357 

[3] AxTELL R., AxELROD R., EpSTEiN J. and Cohen M. D., Computational and Mathematical 

Organization Theory, 1 (1996) 123 
[4] CASTELLANO C, Marsili M. and Vespignani a., Phys. Rev. Lett, 85 (2000) 3536 
[5] Sanchez A. D., Lopez J. M. and Rodriguez M. A., Phys. Rev. Lett, 88 (2002) 048701 



6 



EUROPHYSICS LETTERS 



[6] HoLYST J. A., Kacperski K. and Schweitzer F., Annual Reviews of Computational Physics, 

Vol. IX (World Scientific, Singapore) 2001 
[7] Stauffer D., JASSS, 5 (2002) 1 

[8] Liggett T. M., Interacting Particle Systems (Springer, New York) 1985 

[9] Vazquez F., Krapivsky P. L. and Redner S., cond-mat/0209445, (2002) 
[10] In d = 1 the voter model coincides with the Ising model with Metropolis dynamics at T = 0. In 
higher dimensions it differs from it because of the absence of surface tension; DORNIC I., Chate 
H., Chave J. and Hinrichsen H., Phys. Rev. Lett, 87 (2001) 045701 
[11] Marro, J. and DiCKMAN. R., Nonequilibrium Phase Transition in Lattice Models (Cambridge 

University Press, Cambridge) 1999 
[12] Strogatz S. H., Nature, 410 (2001) 268 

[13] Albert R. and Barabasi A.-L., Rev. Mod. Phys., 74 (2002) 47 
[14] Watts D. J. and Strogatz S. H., Nature, 393 (1998) 440 

[15] Watts D. J., Small Worlds: The Dynamics of Networks between Order and Randomness 

(Princeton University Press, Princeton) 1999 
[16] LiLjEROS F., Edling C. R., Amaral L. A. N., Stanley H. E. and Aberg Y., Nature, 411 

(2001) 907 

[17] Redner S., Eur. Phys. J. B, 4 (1998) 131 

[18] Newman M. E. J., J. Stat. Phys., 101 (2000) 819 

[19] Bartiielemy M. and Amaral L. A. N., Phys. Rev. Lett, 82 (1999) 3180 
[20] Barrat a. and Weigt M., Eur. Phys. J. B, 13 (2000) 547 

[21] Frachebourg L. and Krapivsky P. L., Phys. Rev. E, 53 (1996) 3009; Redner S., A Guide 

to First-Passage Processes (Cambridge University Press, Cambridge) 2001 
[22] AXELROD R., J. of Conflict Resolut, 41 (1997) 203 

[23] Callaway D. S., Newman M. E. J., Strogatz S. H. and Watts D. J., Phys. Rev. Lett, 85 
(2000) 5468 

[24] Cohen R., Erez K., ben-Avraham D. and Havlin S., Phys. Rev. Lett, 85 (2000) 4626 

[25] Pastor-Satorras R. and Vespignani A., Phys. Rev. Lett, 86 (2001) 3200 
[26] D. BOYER and O. Miramontes, Phys. Rev. E, 67 (2003) 035102 



